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We consider the thermally activated escape of an over- 
damped Brownian particle over a potential barrier in the pres- 
ence of periodic driving. A time-dependent path-integral for- 
malism is developed which allows us to derive asymptotically 
exact weak-noise expressions for both the instantaneous and 
the time-averaged escape rate. Our results comprise a con- 
ceptionally new, systematic treatment of the rate prefactor 
multiplying the exponentially leading Arrhenius factor. More- 
over, an estimate for the deviations at finite noise-strengths 
is provided and a supersymmetry-type property of the time 
averaged escape rate is verified. For piecewise parabolic po- 
tentials, the rate-expression can be evaluated in closed ana- 
lytical form, while in more general cases, as exemplified by a 
cubic potential, an action-integral remains to be minimized 
numerically. Our comparison with very accurate numerical 
results demonstrates an excellent agreement with the theo- 
retical predictions over a wide range of driving strengths and 
driving frequencies. 

PACS numbers: 05.40.-a, 82.20.Mj, 82.20.Pm 



I. INTRODUCTION 

The thermally activated escape over a potential bar- 
rier is a recurrent theme in a large variety of physical, 
chemical, and biological contexts In the case of 

foremost practical relevance, the characteristic strength 
of the thermal noise (the thermal energy ksT) is much 
smaller than the potential barrier. As a consequence, 
successful barrier crossings constitute rare events and the 
escape statistics verifies vifith very high accuracy an ex- 
ponential decay as a function of time. In other words, a 
meaningful escape-rate can be defined which completely 
characterizes the decay process. A seminal contribu- 
tion to the theory of escape-rates represents the work 
by Kramers in 1940 which subsequently has been 
refined, modified, and generalized in various important 
directions ||l|-^. 

A particularly challenging novel direction are systems 
far away from thermal equilibrium, either due to non- 
thermal noise or external deterministic forces Q. In 
such a case, the relevant probability distribution strongly 
deviates from the Boltzmann-form in the entire state 
space and its determination becomes a highly non-trivial 
problem. Mutas mutandis, this very same basic diffi- 
culty resurfaces again in all known theoretical methods 
of calculating escape rates in far from equilibrium sys- 
tems H-^. 

The subject of our present paper is one of the sim- 
plest non-equilibrium descendants of Kramers' original 
problem: namely the thermally activated escape of an 



overdamped Brownian particle over a potential barrier 
in the presence of a periodic driving (details are given 
in Sect. II). This is a prototypical setup in the sense 
that investigating the behavior of a system under the 
influence of a periodic forcing represents a particularly 
natural and straightforward experimental situation. Ex- 
amples arise in the context of laser driven semiconductor 
heterostructures [|l^, stochastic resonance [ H , directed 
transport in rocked Brownian motors ]l5|-jl7||, or period- 
ically driven "resonant activation" processes ||l^;0 hke 
AC driven biochemical reactions in protein membranes 
po| , to name only a few. 

Despite its experimental importance, the theory of os- 
cillating barrier crossing in the regime of weak thermal 
noise is still at its beginning. Previous quantitative, 
analytical investigations have been restricted to weak 
|l|-||, slow nil, or fast |2],||,||l driving. In this 
article we continue our recent study p7[ of the most 
challenging intermediate regime of moderately strong and 
moderately fast driving by means of path-integral meth- 
ods. The general framework of this approach is derived 
from scratch in Sect. Ill, thus collecting, streamlining, 
and partially extending previously known material. The 
evaluation of the escape -rate is worked out in Sect. IV 
with the central results ( 116 ) for the time-averaged and 
(108) for the instantaneous escape rate. Especially, these 
results comprise a conceptionally new, systematic treat- 
ment of the rate prefactor multiplying the exponentially 
leading Arrhenius factor. They become asymptotically 
exact for any finite amplitude and period of the driv- 
ing as the noise strength tends to zero. On the other 
hand, for any fixed (small) noise strength, we have to ex- 
clude extremely small driving amplitudes and extremely 
long or short driving periods since this would lead us 
effectively back to an undriven escape problem, which 
is not covered by our present approach. Another situ- 
ation which is excluded in our theory is the case of ex- 
tremely strong driving such that escape events become 
possible even in the absence of the thermal noise |^,^ . 
Closest in spirit to our methodology is the recent work 
p3| , which is restricted, however, to the linear response 
regime (weak driving) for the exponentially leading part 
(Arrhenius factor) and treats the prefactor by means of 
a matching procedure, involving the barrier region only. 
The approximation adopted in that work is complemen- 
tary to ours in that it admits, for a fixed (small) noise 
strength, arbitrarily small driving amplitudes. 

In Sect. V, our analytical predictions are verified for 
the case of sinusoidally rocked metastable potentials 
against very precise numerical results. A first example 
consists in a piecewise parabolic potential, for which our 
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general rate-expressions can be evaluated in closed ana- 
lytical form. In more general cases, exemplified in Sect. V 
by a cubic potential, a few elementary numerically tasks 
remain before actual numbers can be obtained from our 
rate-expressions. The final conclusions are presented in 
Sect. VI. 



II. THE ESCAPE PROBLEM 
A. Model 

We consider the following model for the one- 
dimensional Brownian motion of a particle with coor- 
dinate x{t) and mass m under the influence of a time- 
dependent force field F{x,t): 



mx{t) - F{x{t),t) = -rix{t) + ■ 



(1) 



While the left hand side accounts for the dynamics of the 
isolated particle, the right hand side models the influence 
of its thermal environment with a viscous friction co- 
efficient T] and a randomly fluctuating force ^ [t) , which is 
assumed to be unbiased Gaussian white noise with cor- 
relation 



mm)^5(t-t') 



(2) 



At thermal equilibrium, the intensity D of the noise is 
related to the temperature T according to the Einstein- 
relation D = rjkBT, where fcs is Boltzmann's constant 
101 . Throughout this paper we will restrict ourselves to 
the overdamped motion such that inertia effects mx{t) 
in (|^) are negligible |^^. Choosing the time unit such 
that rj = 1, the stochastic dynamics takes the form 



x{t)=F{x{t),t)+V2Dat) ■ 



(3) 



The force field F{x, t) in (||) is assumed to derive from 
a metastable potential which undergoes an arbitrary pe- 
riodic modulation in time with period T: 



F{x,t + T) ^ F{x,t) 



(4) 



An example is a metastable static potential V{x) as car- 
tooned in fig.l, supplemented by an additive sinusoidal 
driving 



F{x,t) 



-V'{x) 
2 7r/r 



A sin(ri t) 



(5) 
(6) 



Our next assumption is that the deterministic dynam- 
ics in (^) with D — Q exhibits a stable periodic orbit Xs (t) 
and an unstable periodic orbit Xu{t) |3l|, satisfying 



Xs,uit + T) 



F{Xs,u(t)^t) 

■ -^S^U (^) 1 



(7) 
(8) 



where 's,m' means that the index may be either 's' or 
'w'. Moreover, every deterministic trajectory is assumed 



to approach in the long-time limit either the attractor 
Xs{t) or to diverge towards a; = oo, except if it starts 
exactly at the separatrix Xu (t) between those two basins 
of attraction. In other words, the metastable potential 
is required not to be rocked too violently such that par- 
ticles cannot escape deterministically, i.e. without the 
assistance of the random fluctuations in (^. It is clear 
that Xs{t) and Xu{t) must be disjoint and by assuming 
a second "attractor" at a; = cx) we have, without loss of 
generality, implicitly restricted ourselves to case that 



Xuit) > Xs{t) 



(9) 



for all t. Note that the above requirements do not nec- 
essarily exclude the possibility that for certain times t 
the "instantaneous potential" , from which the force field 
F{x, t) derives, does no longer exhibit a potential barrier. 



B. Escape rates 

Next, we return to the stochastic dynamics (|^) with a 
finite but very small noise-strength D such that a parti- 
cle x{t) is able to leave the domain of attraction of the 
stable periodic orbit Xg (t) and subsequently disappear to- 
wards a; = oo but the typical waiting-time before such an 
event occurs is much longer than all characteristic time 
scales of the deterministic dynamics (separation of time 
scales pl p^js^ ). For a quantitative characterization of 
such escape events, our starting point is the probability 
distribution p{x, t) of particles which is governed by the 
Fokker-Planck equation 

Once p{x, t) is known, the population PA{t) of the time- 
dependent basin of attraction A{t) :— {—oo, Xu(t)] of the 
stable periodic orbit Xs{t) follows as 



p{x, t) dx 



(11) 



A suggestive definition of the "instantaneous rate" T(t) is 
then provided by the relative decrease of this population 
per time unit 



m -PA[t)/PA{t) . 



(12) 



We note that particles which leave the domain of attrac- 
tion give rise to a positive contribution to Tit). There 
is also a certain probability that particles from outside 
this domain recross the separatrix Xu (t) , giving rise to a 
negative contribution to r{t). In other words, eq.([l^) is 
the net flux of particles (outgoing flux minus back-flux) 
across the separatrix Xu (t) in units of the remaining pop- 
ulation Pj^{t). By exploiting the deterministic dynamics 
(0) for Xu{t), the Fokker-Planck equation ( |lO| ) forp(x,t). 
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and the definition (O) of PA{t) we can rewrite the in- 
stantaneous rate (|12|) as 



m = 



D dp{xu{t),t) 



PA{t) dxu{t) 



(13) 



Inside the metastable state x < Xu(t) the particle dis- 
tribution is governed by intrawell relaxation processes. 
For small noise-strengths their characteristic time 
scales are well separated from the typical escape time 
itself 1^,^,^. On this time scale of the intrawell re- 
laxation, transients die out and the distribution p(x, t) 
approaches a quasi-periodic dependence on time t. More 
precisely, p{x,t)/ J^"^^^ p{x,t) dx tends, at least for x < 
Xu(t), towards a time-periodic function as t grows. The 
same carries over to the escape probability (^^ and thus 
the time-averaged escape rate 



t+T 



r(t') dt' 



(14) 



becomes independent of the time t. 

Our assumption of weak noise guarantees that the loss 
of population fU(i) is negligible on the time scale of the 
intrawell relaxation for any initial distribution p{x,to) 
that is negligibly small in the vicinity and beyond the in- 
stantaneous separatrix a;„(to). The denominator in dl^ ) 
can thus be approximated by 1 for all times t — to much 
smaller than the characteristic escape time 1/f itself. 
Further, we can restrict ourselves to delta-distributed ini- 
tial conditions of the form p(x,to) ~ S{x — xq) with xq 
inside the basin of attraction A{t) of Xs {t) such that the 
overwhelming majority of realizations (0) will first re- 
lax towards a close neighborhood of the attractor Xs{t) 
before they escape. The behavior of more general ini- 
tial distributions then readily follows by way of linear 
superposition. Moreover, one expects ||l|,^,p2j that af- 
ter transients (intrawell relaxation proccsses)liave died 
out, the time dependent escape rate ( p^ ) will actually 
become independent of the initial conditions xq and Iq. 
Denoting by p{x, t\ xq, io) the conditional probability as- 
sociated with an initial delta peak at xq we thus can 
rewrite (|3|) as 



m = -D 



dp{xu{t),t\xa,to) 
dxu(t) 



(15) 



We recall that this expression is valid even ii t — to is not 
large, but then T{t) still depends on xq and Iq. On the 
other hand, t — to has been assumed to be much smaller 
than the typical escape time l/T. However, on this time- 
scale the rate T{t) has practically converged to its asymp- 
totically periodic behavior and thus the extrapolation of 
r{t) to arbitrarily large t — to is trivial. 



C. Supersymmetry 

Given a time-periodic force field F{x,t), we define its 
supersymmetric partner field F{x, t) [ ^2|j3^j35[] according 
to 



F{x,t) -.^Fi-x^-t) 



(16) 



For instance, if the force field F{x,t) derives from a pe- 
riodically rocked potential V{x) like in (^), then its su- 
persymmetric partner is obtained by turning V{x) up- 
side down, followed by an inversion of the x-axis, i.e. 
V{x) = —V{—x), see fig. 2, while the driving A sm{ilt) 
in (^) remains invariant (up to an irrelevant phase). In 
such a supersymmetric partner field F{x,t), the stable 
and unstable periodic orbits exchange their roles, thus 
defining a new escape problem out of the basin of attrac- 
tion of the new stable orbit Xs{t) = —Xu{—t) across the 
new separatrix a;„(i) = —Xs{—t). 

It has been demonstrated in |£2|,|^ that for force fields 
F{x,t) of the form (||), the time-averaged rate ( |l^ is in- 
variant under the supersymmetry transformation ([l^). 
The same line of reasoning can be readily gener- 

alized to force fields of the form F{x, t) ~ —V'{x) + y{t) 
with an arbitrary periodic driving y{t). In our present 
paper we will show that for asymptotically weak noise D 
the time averaged escape rate ( |7^ is invariant under the 
general supersymmetry transformation ( [7^ without any 
further restrictions on F{x,t). Regarding the notion of 
supersymmetry and especially its connection with super- 
symmetric quantum mechanics, we refer to p^js^ and 
further references therein. We finally remark, that the 
standard definition of the supersymmetric partner force 
field is —F{x, —t). For our present purposes, the defini- 
tion ( p^ ) is equivalent but more convenient. 



III. PATH-INTEGRALS: GENERAL 
FRAMEWORK 

In this section the general framework of a path-integral 
approach to the stochastic dynamic s ^) is outlined. 
Though these concepts are not new [p|- |12| , |36| -^ , we find 
it worth while to briefly review them here in order to 
make our paper self-contained. We also note that most 
of this section remains valid beyond the particular as- 
sumptions on the force field F{x,t) from Sect. II. 



A. Time-discretized path-integrals 

Much like in quantum mechanics, also in the present 
context of stochastic processes, path-integral concepts 
have a tangible meaning only when considered as the lim- 
iting case of appropriate discrete-time approximations. 
Our first step is therefore a discretization in time of the 
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overdamped stochastic dynamics (||). Denoting the ini- 
tial and final times by to and tf, we introduce the defini- 
tions 



At 



= [tf - to]/N 
= to + nAt 



(17) 
(18) 
(19) 



where n = 0,1, ...,7V. The integer N is considered as 
large but finite and will ultimately be sent to infinity 
(continuous time limit). The discretized dynamics (|^) 
then takes the form 



F{Xn,tn)At + V2DAt^„ 



(20) 



where the ^„ are independent, identically distributed 
Gaussian random numbers with probability distribution 



,-1/2 



(21) 



As a side remark we notice that the so-called "prepoint 
discretization scheme" pO|-|4^ (not to be confused with 
the Ito-scheme in the stochastic dynamics (||)) has been 
implicitly adopted in ( pO[ ) for the sake of later conve- 
nience. Other "discretization schemes" |40[ - |42[ would 
give rise to a somewhat modified path-integral formal- 
ism but would - of course - lead to identical results as 
far as the actual stochastic dynamics (^ is concerned. 
In passing we further note that our treatment for Eq. (|^) 
can be generalized to multiplicative noise g{x) C{t)i with 
g{x) 0, without encountering additional difficulties. 

For the conditional probability pNixn+i,tn+i |a;„,i„) 
to reach the point Xn+i at time tn+i when starting out 
from Xn at the previous time step i„ we find from the 
discretized dynamics ( pO| ) and the noise-distribution ( |2l| ) 
that 

PN{x,i+l,tn+l I Xn, tn) = 

J S{Xn+l - Xn - F{Xn,tn)At - V2DAt^n) P{&i) d^n 
1 ( [Xn+1 - Xn~ F(x„, t„) A<]2 



{4:7: D At) 2 



exp 



ADAt 



where 
Sn{xo, ■■■,xn) 



Here and in the following, integrals over the entire real 
axis are written without the integration limits ±oo. Fur- 
ther, the mutual independence of the random numbers 
^„ in (^0|) (Markov property) implies for the conditional 
probability the Chapman-Kolmogorov relation 

PN{Xn+2,tn+2 \ Xn, tn) = 

PNiXn+2,tn+2\Xn+l,tn+l) PN{Xn+l,tn+l \ Xn,tn) dXn+1 



Upon iteration of this relation in combination with 
( p^ ) one finds for the conditional probability the time- 
discretized path-integral representation 



PN{xf,tf I Xo,to) = 

dxi ■ ■ ■ dxN- 



I 



(4 7rDAi)^/2 



exp 



SNixo, 



,xn) 



D 



(23) 



2-^ A 



At 



(24) 



is the discrete-time "action" or "Onsager-Machlup func- 
tional". While xi, x^r-i are integration variables in 
(p3|), the initial- and end-points are fixed by the pre- 
scribed a;o and by the additional constraint xjv ~ Xf, see 



B. Saddle-point approximation 

For small noise-strengths D the path-integral (|2^) is 
dominated by the minima of the action Sn^xq, ...,xn)- 
The existence of at least one (global) minimum can be 
readily inferred from the general structure of the action in 
(p^. To keep things simple we assume for the moment 
that besides this global minimum no additional (local) 
minima play a role in (p3|). Denoting the global mini- 
mum by X* := {xq, xf^ it follows that it satisfies the 
extremality conditions 



dSN{^*) 
dx* 



= 



(25) 



for ?T.= 1,...,A^— 1, supplemented by the boundary con- 
ditions for n = 0,N: 



Xo 



(26) 



Under the assumption that the noise-strength D is small, 
the path-integral in ( p3|) can be evaluated by means of a 
saddle point approximation about the minimizing path 
X* with the result 



(22) 



PN{xf,tf I Xo,to) = ZAr(x*) e 



-Sjv(x*)/D 



(27) 



where the prefactor Z7v(x*) is given by a Gaussian inte- 
gral of the form 



ZAr(x*): = 



dyi--- dyN-i 
(4 7rDAi)^/2 



X exp ■ 



2D ^ ^" 



n.m — 1 



d^Sjx*) 
dx*ndx*^ 



(28) 



and where in the order of magnitude expression 0{D) 
only the dependence on the noise-strength D is being 
kept. A more detailed quantitative estimate of this cor- 
rection 0{D) is a difficult, and to our knowledge unsolved 
task. 

The Gaussian integral in ( |2^ ) is readily evaluated to 
yield 
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Zn{^*) 



AnDAt det 2 At 



(29) 



where det(A„m) indicates the determinant of an iV — 1 x 
TV — 1 matrix with elements Anm- As demonstrated in 
Appendix A, the determinant appearing in ( |2^ ) can be 
rewritten in the form of a two-step (second order) Hnear 
recursion 

Qn+l ~ 2 Q* — Qn-l 



At^ 

01 i„) - q;_i ^^'(x;_i, t„_i) 



At 



At 

F{x^, tn) 



with initial conditions 

g* = Ai , 9l_^ = i + oiAt) 
from which the prefactor Zn{x*) in (E9) follows as 



(30) 



(31) 



1-1/2 



(32) 



The fact that x* is a minimum of the action ( p4| ) guar- 
antees that > 0. Here and in the following we use 
the abbreviations 



F'ix,t) 



dF{x,t) 
dx 



■ , , dF(x,t) , , 
F(x,t):=^^ (33) 



and bracket-saving expressions like /(x)^ are understood 
as {!{x)f. 

As we shall see later, we have to leave room for the pos- 
sibility that even for small noise-strengths D more than 
one (global or local) minimum of the actio n (|2^ ) notably 
contributes to the path- integral expression (23). We label 
those various non-negligible minima by the discrete in- 
dex k but leave for the moment the precise set of fc-values 
unspecified. Each of the minimizing paths xl thus sat- 
isfies an extremality condition of the form (pS]). Under 
the assumption that those minima x^ are well separated 
in the N — 1-dimensional space of all paths (xq, ...jXat) 
appearing in (|2^), the saddle point approximation ( |2^ ) 
simply acquires an extra sum over fc with a correspond- 
ing extra index k in Combining (|2|) and (||) 
we thus arrive at 



PN{xf,tf I XQ,to) = ^ 



-S«(x*)/D 



r [1 + 0{D)] . (34) 



C. Continuous-time limit 

Next we turn to the continuous-time limit N ^ oo, 
Ai — > in (^7|). The continuous-time conditional prob- 
ability p{xf,tf I xo,to) when TV ^ oo in ( p3| ) is symboli- 
cally indicated by the path-integral expression p& 



x{tf)=Xf 

p{xf,tf\xo,to)= J I?x(t) e-^[-(*)l/^ 

3;(to)=a;o 

where 



(35) 



S[x{t)] 



L{x(t),x{t),t)dt 



(36) 



to 



is the continuous-time limit of the action (p3) with 



L{x, i, t) 



F{x,t)f 



(37) 



as Lagrangian. The extremality conditions for the min- 
imizing paths x'^{t) in the continuous-time limit are ob- 
tained from (U) and (|2|) by letting At — > as 

= F{<{t),t) + F{xl{t), t) F'ixlit), t) (38) 

with boundary conditions (cf. (p6|)) 

xlito) = , xl{tf) = Xf . (39) 



The same result (^8|) can also be recovered as the Eulcr- 
Lagrange equation corresponding to the Lagrangian ( |3^ ) . 

Equivalent to this Lagrangian dynamics is the follow- 
ing Hamiltonian counterpart: 



H{x,p, t) :— px — L = p^ + pF{x, t) 

m) = -Pi{t)F'{xi{t),t) 

xUt) = 2pl{t) + F{xl{t),t) . 



(40) 
(41) 
(42) 



The last equation (|4^) may also be considered as the defi- 
nition of the momentum {t) in terms of {t) and {t) . 
With (|3^) and (^) the action (^6|) along a minimizing 
path a;^(t) follows as 



</>fc(x/,t/) ■.^S[xl{t)]= f pHtfdt 

J to 



(43) 



where the dependence of the action (j)k{xf,tf) on the ini- 
tial condition xq at time tg has been dropped. For later 
use we also recall the well-known result from classical 
mechanics that the derivative of the extremizing action 
with respect to its endpoint equals the canonical conju- 
gate momentum, i.e., 



bk{Xf,tf) 

dxf 



pUtf) . 



(44) 



Finally, the continuous-time limit for the conditional 
probability (B4h in combination with (Esh takes the form 



p{xf,tf\xo,to) = 



^-<t>kixf,tf)/D 

[4nDQl{tf)]-2 



[1 + 0{D)] , (45) 



where Q^{t) is governed by the second order homoge- 
neous linear differential equation l27|j39| , ^ that follows 
in the limit At ^ from ^ and (^JT 
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+ QUt)pl{t)F"{xl{t),t)^0 



(46) 



Similarly, the initial conditions ( pi] ) go over for At ^ 
into 



Qlito) = , QUto) = 1 . 



(47) 



We remark that according to (|38| ) and (|3^) the mini- 
mizing paths x'l (t) are independent of the noise-strength 
D. Consequently, neither 4>k{xf,tf) from ( ^ , p9 43) nor 



I.e., 



Ql{t) from (||,|^ depend on the noise-strength 
no implicit additional /^-dependences are hidden in (45) 
We further note that by means of the substitution 



In view of ([45[), the instantaneous rate ( |5l| ) has the sug- 
gestive structure of "probability at the separatrix times 
velocity" . 

As already mentioned in Sect. II, for sufhciently large 
times tf — to the instantaneous rate ( |5l| ) is expected to 
become independent of the initial position xq as long as 
xq is located inside the domain of attraction of the stable 
periodic orbit Xs{t). A more detailed discussion of this 
point will be given in Subsect. IV. G. To keep things as 
simple as possible we focus in the following subsections 
on the particular case that Xq is located at the stable 
periodic orbit, i.e., 



Xo = Xsito) . 



(52) 



9l{t) 



2Ql{t) 



F'ixUtU) 



(48) 



the linear homogeneous second order equation ( pq ) goes 
over into the non-linear first order Riccati-equation 



gl{t) + 2gUt)^ + 2gl{t)F'ixl{t),t) 
-pUt)F"{xUt),t) . 



(49) 



Since (|4^) does not lead to a meaningful initial condition 
for gl{t) in (^), the Riccati-equation ( |49| ) can only be 
used for times t > t^. For this reason and also from the 
viewpoint of calculational efhciency we found that for 
practical purposes the linear second order equation ( ^ ) 
is often superior to the Riccati-equation (p9|). 

To establish contact with previously known results we 
finally remark that one can identify 



glitf) 



(50) 



This relation (|5^) and the associated Riccati-equation 
(^) are usually derived by introducing a WKB-type 
ansatz into the Fokker-Planck-equation for the condi- 
tional probability distribution (cf. (|l^)) and then com- 
paring powers of the noise-strength D. Since a direct 
derivation by means of path-integral methods is not 
known to us, we have included such a derivation of the 
relation (50) in the Appendix B. 



PATH-INTEGRAL SOLUTION OF THE 
ESCAPE PROBLEM 



IV. 



By introducing the path-integral expression (^5|) for 
the conditional probability into the formula (15) for the 
instantaneous rate Tit) at time t = tf and taking into 
account (El) we obtain our first main result p7|] , namely 



ntf) = E 



[A^DQl[tf)Y/^ 



[1 + 0{D)]. (51) 



A. Minimizing patiis 

Our next goal is the characterization of all the mini- 
mizing paths x*f^{t) which significantly contribute to the 
sum in (51). Our first observation is that for any finite 
to and tf the action (36) exhibits in the generic case a 



unique global minimum respecting the boundary condi- 
tions 



4(io) ^ Xs{to) , Xl{tf) = Xuitf) 



(53) 



according to ( |39| ) and ( piy52| ). To be specific, we denote 
this globally minimizing path as [t) . From the explicit 
form of the Lagrangian (37) we can infer that for large 
values oi tf — tQ the minimal path x*f.^ (t) follows most 
of the time rather closely a deterministic trajectory, i.e., 
(t) ~ F{xl {t),t), in order not to accumulate a too 



large amount of action (p^). In view of (|^) and (53) it is 
thus suggestive that xl^ (t) starts at xl^ {to) = Xsijo) and 
then continues to closely follow the stable periodic orbit 
Xs{t) for quite some time. At a certain moment, x^^it) 
leaves this neighborhood and travels in a comparatively 
short time into the vicinity of the unstable periodic orbit 
Xu(t), where it remains for the rest of its time and ends 
at x*f.^{tf) = Xu{tf). Only during the crossover from the 
neighborhood of Xs{t) into that of Xu{t) does the path 
x^^ [t) substantially deviate from a deterministic behav- 
ior and so gives rise to the main contribution to the ac- 
tion (p6[). We desist from a more rigorous derivation of 
these basic qualitative features since they are quite simi- 
lar to the well-known barrier-crossing problem^n a static 
potential (time-independent force-field) 
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Espe- 
cially, the relatively short "crossover-segment" of x\^ {t) 
between the long sojourns close to the stable and unsta- 
ble orbits has lead to the name "instanton" for such a 
path. 

As we will see in more detail later, a meaningful limit 
of x^.^ (t) exists for to —oo and — > cxd (henceforth ab- 
breviated astf — to — > oo) in the sense that xl^{t) follows 
closer and closer the periodic orbits Xs,uit) over longer 
and longer time intervals, while the "crossover-segment" 
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does practically not change its shape any more. Also the 
associated minimal action S[xl^{t)] from (|3^) tends to a 
finite limit. In fact, one can readily show that the mini- 
mal action cannot increase upon increasing tf and/or de- 
creasing Iq. Since it is furthermore bounded from below, 
the existence of the limit follows for the action as well as 
for the minimizing path itself. More iinportantly, from 
the time-periodicity of the force field (||) one can infer 
that in the limit t f — to ^ oo the action S[x'l^{t + nT)] 
has the same value for any integer n. In other words, 
for infinitely large tf — the action no longer exhibits a 
unique absolute minimum, rather each path x'^^{t + nT) 
globally minimizes the action. However, these degener- 
ate absolute minima are still well separated in the space 
of all paths x{t) appearing in (l35|). This feature is the 
salient difference between our present problem and its 
time-independent counterpart ||37| , ^ , ^4| -|4^ , which ex- 
hibits a continuous degeneracy (Goldstone mode) in the 
limit tf — to —^ oo. Put differently, the time-periodic 
force field reduces the continuous time-translation sym- 
metry into a discrete one. Since the rate-formula ( |5l| ) 
assumes well separated minima a;^(i) of the action, it is 
quite clear that the time-independent case must be ex- 
cluded in the following. 

We emphasize that the minimizing paths xt{t) remain 
well separated and thus the rate-formula becomes 
asymptotically exact for any (arbitrary but fixed) finite 
values of the driving amplitude and period as the noise 
strength D tends to zero. Apart from this fact that 
in the limit D ^ the 0{D) correction in the saddle 
point approximation ( p7| ) and thus in (^) vanishes, a 
more detailed quantitative statement seems difficult. On 
the other hand, for a given (small) noise strength D, 
we have to exclude extremely small driving amplitudes 
and extremely long or short driving periods since this 
would lead us effectively back to the static (undriven) 
escape problem, which requires a completely different 
treatment (especially of the (quasi-) Goldstone mode 
|||,|^J38[||-|§) than in (|3). Put differently, in any 
of these three asymptotic regimes, the error 0{D) from 
(^,^ becomes very large. 

For later reference we denote the minimizing path 
^fco(^) when tf — to oo by x*p^{t), keeping in mind 
that we are still free to shift its time argument by an 
arbitrary multiple of T. The corresponding action is 



fopt 



= S[x*t{t)] = ^ lim 



mm 

x(t) 



S[x{t)] , (54) 



where the second identity may also be considered as an 
implicit definition of x*p^{t). Similarly, any other quan- 
tity associated with x*p^{t) will be marked by an index 
"opt", for instance p*pt(t) (see (g)), Ql^.it) (see (|6|)), 
andg*pt(t) (see (|8|)). 

In principle, besides the absolute minimum x*p^{t) of 
the action there may coexist further (absolute or relative) 
minima which cannot be identified with each other after a 



time-shift by an appropriate multiple of T. While the co- 
existence of further absolute minima is non-generic, coex- 
isting relative minima are irrelevant for sufficiently small 
noise-strengths D as far as the sum in (|5^) is concerned. 
Though both cases could be easily taken into account in 
the following discussion, we will restrict ourselves to the 
simplest and most common case that x* f it -\- nT) are 
the only (relevant) minima of the action (M) in the limit 



ti 



to 



Returning to finite but large values oftf — to, we expect 
- as a precursor of the tf—to oo limit - that besides the 
unique absolute minimum x^^ (t) there will coexist many 
additional relative minima x^(t) with an only slightly 
larger action. All those minima xl{t) possess a limit 
when tf—to ~^ oo in the same sense as for the case k = ko 
described above (quantitative details will be given later). 
Moreover, when tf — to ^ oo then each x\{t) approaches 
x*p^(i H-ri(fc) T) for a suitable choice of n{k) and without 
loss of generality we can assume a (re-)labeling of the 
x*i.{t) such that n{k) — k. In other words, to each x\{t) 
belongs a very similarly looking ''^master-path!^ ^opti^ + 
kT), see fig. 3. Since tf—to is finite, there is a finite 
number (of the order (t/ — to)IT) of minimizing paths 
a;^(t) and without loss of generality we can assume that 
the indices in the sum (O) start at A; = and run until 
a certain maximal value K{tf,to)'. 



0<k<K{tfM)^0{{tf-to)/T) 



(55) 



Thus Xo{t) is that minimizing path which closely follows 
Xs [t) as long as possible and crosses over to the neighbor- 
hood oixu{t) at "the latest possible moment" (see fig. 3), 
and similarly for 2^x(tj: to)^^^' 

Note that all the general qualitative features discussed 
above are nicely illustrated by the explicit example in 
Sect. IV.A. 



B. Neighborhood of periodic orbits 

Our final goal is to approximate the action (f)k{xf,tf) 
and the prefactor pl.it f)/[Ql.it f)]^^^ for all minimizing 
paths xl.{t) that play a non- negligible role in the sum 
( ^ , p5| ) solely in terms of the "master-path" a;*pj(t) and 
its descendents plpt{t), Qoptl^): ^tc. To this end we first 
address the behavior of a path x^t) within the neigh- 
borhood of either the stable periodic orbit Xs{t) or of 
the unstable one Xu{t). Within these regions, the time 
dependent force field can be approximately written as 

F{x,t) = Fix,,uit),t) (x - xs,uit)) F'{xsAt)^t) ■ 

(56) 

Note that these approximations are valid not only if 
X — Xs,u(t) is small but also if F"{y,t) is small for all 
y between x and Xs^„(f). An immediate consequence of 
(p6|) are the relations 
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F'{x,t)=F'{xs,u{t),t) , F"{x,t)=0 



(57) and similarly 



As long as a minimizing path remains in a re- 

gion where these approximations apply, the Hamiltonian 
equations (41 14^) take the form 

(58) 
(59) 

(60) 

(61) 
(62) 



AxUt) = 2pUt) + Axlit) F\xsAi), t) , 
where we have introduced 

Axlit) := xlit) - Xs,uit) ■ 
Their solutions are: 

Axlit) = A4(ii)eA-"(*'*^) +plit)Is,uit,ti) 



where ti is an arbitrary reference time (within our as- 
sumption that (j5^) applies for all the considered times t) 
and where 



^s.uit, ti 



F'ixs,uit').t')dt' 



'ti 



(63) 
(64) 



Obvious properties of the functions As^uit,ti) from 
( |63| ) are: 

A,,„(t,i2) = A,,„(i,ti) + A,,„(ti,t2) (65) 

A.,„(ti,t) - -A,,„(t,ii) (66) 

A,,„(t + T,ti+T)^ A,,„(t, ti) . (67) 



Further, one readily sees that the quantities 
A.,„ := A,^„(t + r,t)/T 



(68) 



are indeed <-independent and characterize the stabil- 
ity/instability ("Lyapunov exponents") of the periodic 
orbits, namely 



A., < 



A„ >0 



(69) 



One even expects that Asit, ti) < and A„(t, ti) > not 
only for t-ti = T,2T, ... (cf .^)) but in fact for all t - 
ti > 0, however exceptions cannot be excluded for not too 
large t — ti. From ( |6^ ) and ( |68| ) it follows that As,„(t,ti) 
can be written as the sum of a linear function X^^u 'it~ti) 
and a periodic function of t. As a consequence, we obtain 



As,ti(i,ii) ~ Xs,u ■ it — t]) 



(70) 



for asymptotically large positive and negative times t~t\. 

Turning to the discussion of Is,uit,ti) from (^, we 
first note that 



Isit,t,) = /,(*)- e2A=(*^*i)/,(ii) 



Isit):^ lim Isit,ti) 



(71) 

g2A.(t,t')rfi' (72) 



Iuit,ti) = -/„(t) + e2^"(*'*^)/„(ii) 



(73) 



luit) := - lim Ut,ti) = 2 / e2^"(*'*') dt' . (74) 

ti — >oc 



Thus, Is,uit) are positive and finite for all t and obey 



/.,«(t + T) 



(75) 



It follows that Isit, ti) in (|7l|) is given by a periodic func- 
tion of t minus the product of another periodic function 
times an exponentially decreasing factor exp{As • (t— ii)}, 
and analogously for /„(i,ti) in (|73|). 

Choosing as reference time ti — to in ( |6^ and tak- 
ing into account that AxUto) = according to ( p2||60| ), 
implies that in the neighborhood of Xg it) we have 



Axlit) =plit)Isit,to) 



(76) 



Dividing this result by the same identity evaluated at a 
different reference time ts > to and taking into account 
( |6l| ) we obtain 



A4(t) = A4(^.)e-^=(*-*=)^^^ (77) 

(78) 



p^(t) = A4(i,)e-^=(*-*=V^.(i.,to) 



Both these expressions consist of an exponentially in- 
creasing factor exp{— Ag ■ it — t^)} times some periodic 
function of t. In ( |77| ) one has in addition a quickly de- 
creasing correction. The corresponding behavior in the 
neighborhood of x„(i) is given by 

Axlit) ^ Plit) hit,tf) (79) 

Axlit) = Axlit^) e-^"(*-*") (80) 



Pk 



(t) = A4(t„)e-^"(*'*")//„(t„,i/) , 



(81) 



where t„ is some reference time with t^ < tf. As ex- 
pected, (|8^,|l|) are now dominated by an exponentially 
decreasing behavior exp{— ■ it — t^)}. We further re- 
mark that for the master-path x*„|.(t) we have to —00 
and tf ^ 00, thus Is uit^o.f) in ( |76[]8l| ) go over into 
Is.njt) according to ( |7l|j7^ ) and so all four equations 
(0,^,|o|j8^) are exactly given by exp{-As,„ • (i - t^.n)} 
times certain periodic functions of t. 

Within our above local analysis of the neighborhoods 
of Xs,uit)^ the reference times tg^u are still arbitrary and 
the corresponding parameters AxUtg^u) remain undeter- 
mined. They can only be fixed through the global behav- 
ior of xlit). It is instructive to reconsider the same thing 
from a somewhat different viewpoint. From ( [7l| , [72|j7^ ) 
we conclude that 



PliU) 



to 

= hits) -2 J 



=2 A. di' 



(82) 
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and similarly 



Plitu) 



(83) 



starting point t^. Formally, this approximation is equiv- 
alent to letting 



<o — oo 



(85) 



Let us consider ts > to as fixed and such that the approx- 
imation (^6|) is valid for all t G \to,ts\- Within the same 
restriction, we now consider the quantity Ax^(ts) as a 
parameter. For any value of Ax*^{ts), eq.(^) thus fixes 
pl{ts)- With these initial conditions for x^t) and pl{t) 
at time t = ts on may then propagate the Hamiltonian 
equations ( 4l]j4^ ) up to the time t = tf. It is clear that for 



a typical c 



loice of Axl{ts) such a "shooting procedure" for t — > -cx) and a;*pt(i) 



C. Approximations in terms of the master path 

Our next objective is to express the action ( p3| ) of the 
path xl{t) in terms of the associated master path xlp^{t+ 
kT). We recall that while a;^(i) satisfies the boundary 

Xs{t) 

We now 



conditions (p3|), those of a;*pj(i) are a;*pj(t) 

Xu{t) — > for t 



oo 



does not lead to the desired end-result Axl{tf) = 0. But 
we also know from the mere existence of the minimizing 
paths that there must be specific Aa;^(ts)-values which 
do the job. Furthermore, eq.(83) tells us that it is not 



necessary to proceed until t — tf, rather it is sufficient to 
take any time i„ at which xl.{t) has reached the Xu{t)- 
neighborhood and then check whether (|3|) is satisfied. 

If ts ~ to and tf ~ tu are already large then the inte- 
grands in (|2|,|3|) are extremely small. Thus, tiny changes 
of Axl{ts) and p%{ts) will lead to huge changes of t^ 
and tf. Especially, by letting t^, —oo and tf oo 
those integrals vanish and one recovers the master path 
x*p^(a; -|- kT) associated with x*^,{t). This confirms our 
conclusion from the previous subsection that a meaning- 
ful limit of {t) for t^ ^ oo and tf ^ oo exists and that 
for finite but large tf — to the difference between a;^(i) 
and the associated master path x*p^ {t + kT) is extremely 
small for alH S [io,i/]- 

An example for which tf ~ t^ is not large is the path 
Xo{t), i.e., the one which crosses over from the neighbor- 
hood of Xs (t) into that of x^ (t) at the latest possible mo- 
ment, see fig. 3. For this path Xo(t), the time t„ at which 
it enters the neighborhood of Xu (t) is already rather close 
to tf an so the integral in ( ^ ) is not any more small. As 
a consequence, the deviation of Xo{t) from x*pj(t) is no 
longer small as t approaches tf. In particular, for t = tf 



it follows that Xn 



''opt 



(tf) 



-^Kptitf) is no longer 



small and with (f79) we conclude that the same is true for 
the momentum p*pj(i/), i.e. 



P*opt(tf) 



not small 



(84) 



With increasing fc- values, the deviations —Axlp^{t + 
kT) between xl{t) and the associated master path 
x*pt(t -I- kT) in the vicinity of are rapidly decreas- 
ing, essentially like exp{— A„A:T}, see ( p8pO| ). In the 
same way, for the largest possible fc- values, k ~ K{tf,to) 
(see (11])), corresponding to paths x^t) with only a very 
short initial time segment close to Xs{t), the deviations 
from a;*pj (t -t- fc T) are no longer small for t close to to. As 
we will see later, paths a;^(i) with such large fc- values are 
negligible in the sum ([sl]). For this reason, we will hence- 
forth neglect deviations between x^ (t) and a;*pt (t + kT) 
and between p^, {t) and p*^^ {t + kT) for times t near the 



modify the latter boundary condition and require instead 
that 



tk -.^tf + kT 



(86) 



is the final time and Xk 



^opt 



(tk) the final position. In 



other words, we simply truncate the master path a;*pt(^) 
at the time tk, associated with the final time i/ of xl.(t). 
Since this "new" path x*p^(t) with t € [—oo,tk] obviously 
still satisfies the Hamiltonian equations (^ 42 ) it is again 
an extremizing path. The value of the action for this path 
follows like in ([43) as 



(kaptixk^tk) := / P*opt{t) dt 



and the relations (^,^^ take the form 
d4>optixk,tk) 



dxk 
d'^ (t>opt{xk,tk) 
dxl 



Popt(ifc) 



9*opt{tk) ■ 



With (|3|j54|) we can rewrite (|87|) as 

(t>opt(.Xk,tk) = (/>opt - / Popt(i)^ dt 



(87) 

(88) 
(89) 

(90) 



Next we express the action (36) of the path x*f.{t) by 
expanding the one belonging to the associated master 
path x*^p^[t + kT) in powers of the difference —Ax'^p^{tk) 
between the end-points x*f.{tf) = Xu{tf) and a;*pt(i/ -|- 
kr) ^ xlp,{tk): 



'l^k{xl{tf),tf) = (j)opt{xk:tk) - Ax*{tk) 



d4>opt{xk,tk) 



dxk 



A<pt(ifc)-' (t)opt{xk,tk) 



dxl 



(91) 



As justified above eq.(|85|), the analogous contribution in 
powers of Aa;*pj [to -f kT) is negligible on the right hand 
side of (91). By exploiting (^-90) and the counterparts 
of (|79-81) for a;*pt(t-t- kT), one arrives after a short cal- 
culation at 
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Jtk 

X [I + g^pdtk) Iu{tk) + 



(92) 



A similar expansion of p*pt(tfc) from (g8|) yields for pl{tf) 
the approximation 



Plitf) = Plptitk) + ^xlptitk) 



dxl 



+ ... 



" -Popt 5opt {tk)Iu{tk) + 



(93) 



We now turn to the prefactor-terms QX{t) in (51). 
Within the neighborhoods of Xg u{t) for which the ap- 
proximations (|5^) and thus ( ^7|) are valid, we can infer 
from (E|) that 



Ql{t)/2 ~ QUt) F'{x,,u{t), t) = const. . (94) 

By comparison with (Eq) we further see that 



(95) 



The constant /i^ , which is connected with the neighbor- 
hood of Xs{t) and is typically different form follows 
from the initial conditions ( [47|) as /is = 1/2. Hence, the 
solution of (M) takes the form 



):.(t) = /,(t,to)/2 



(96) 



As a by-product we find from (|50|), evaluated for an ar- 
bitrary final condition tf — t and Xf = x in combination 
with (pSpl) that 



dx"^ 



Isit,to) 



(97) 



Within the linearization (|5^), closer inspection shows 
that only a single summand appears in the conditional 
probability ( p5| ) and one recovers the expected Gaussian 
result for x close to the stable periodic orbit Xg (t) : 



p{x,t\Xs{to),t()) 

1 

2 7rDIs{t,ta) 



cxp 



2DIs{tM) 



(98) 



Returning to (P6|), it is remarkable that besides the 
initial time no further details of the path a:^ {t) play a 
role. Especially, if x*j.{t) remains for a long time in the 
neighborhood of Xs(t) where (^) is valid, then by the 
time it leaves this neighborhood, say t = ts, the quantity 
Is{t,to) is practically equal to hit) from ( [72[ ) and thus 
Ql{t) equal to the associated master-prefactor Qopt(^ + 
k T) . Within our usual approximation (pSf) we thus have 



(99) 
(100) 



These relations are then used as initial conditions in 
(|^) in order to propagate QX{ts) and Q*pt(4 + kT) 



through the "crossover-segments" of the corresponding 
paths xl.{ts) and a;*pt(is + kT) up to a certain time- 
point, say t — tu, beyond which the linearization ( |5^ ) 
about Xu{t) and thus (|9^ ) can be applied. 

Once the neighborhood of Xu{t) is reached, i.e., for 
t > tu, the solution of ( p4|j9^ ) can be written with (jsj) 
as 

Q:.(t)=Q:(t„)e2^"(*^*")[l-.9:(i„)/.(t„,t)] . (101) 
In view of (pW this yields furthermore 



roptit)p:pt{t)' ^ Q*optitu)p*opt{tu)' 

X 9opt{tu)Iu{tu,t)] 



(102) 



Due to ([73|), the factor /„(*„, i) approaches —Iu{tu) as 
i — becomes large. It follows that the left hand side of 
(102) tends towards a finite limit as i ^ oo: 



gopt lim Q*{t)p*{t) 



(103) 



Since i„ is an arbitrary reference- time in (102), we can 
first let t — > oo and then rename <„ as t with the result 



gopt 

piA^y 



fJ-opt Iu{t) 



(104) 



where the (finite) constant /iopt is defined analogously to 
m as 



/iopt := hm 5*pt(<)g*pt(t) 



(105) 



Exploiting once more (^5|) and (105), we can eliminate 
Qopt {t) in ( |104D in favor of 5*pt (t) with the result 



glptit) 



plpM 



qopt/ ^J■opt - Popt{t)'^Iu{t) 



(106) 



As discussed below (|83|), the deviations of a;^(t) from 
the associated master path a::*p^ {t + kT) become smaller 



and smaller as k increases and in view of (|99|, 100|, 46|) we 
expect a similar convergence of Ql{t) towards Qopt(^ + 
kT). In Appendix C, the following quantitative estimate 
for this convergence is established for all times t G [tu,tf]: 



QUt)^Q*dt + kT) [l + 0{p*{tk)')] 



(107) 



where the order of magnitude is meant with respect to 
the dependence on k. 



From t he te chnical viewpoint, Eq.(104) in combina- 
tion with (107) is a central and highly non-trivial result 
of our present work. Since /„ (t) is periodic in t and since 
p*pt (i) decreases exponentially according to (^TJ), we see 
from ( |104 107) that the prefactor Ql.{t) diverges expo- 
nentially with the time which the path x^, (t) spends in the 
neighborhood of Xu{t), in striking contrast to the behav- 
ior ( p6| ) close to Xs{t). The basic physical reason for this 
divergence of Ql{t) is that the probability of a stochastic 
process (0) to permanently remain close to the unstable 
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periodic orbit {t) decreases exponentially with increas- 
ing time. Since typically the process closely follows a de- 
terministic trajectory, the action barely grows and thus it 
is the prefactor 1/Q^(i/)^^^ in ( ^5|) which has to account 
for the exponential decrease in time! 

Since Popt(^fe) decreases exponentially with k we see 
from (106) that 5opt(^fe) tends to zero like p*pt(tfe)^. In 
view of ( |88| , p^ ) we therefore conjecture that also higher 
derivatives of (poptixk^tk) continue to scale like the cor- 
responding powers of p*pj(ifc). The terms indicated by 
the dots in (|9l[-p3|) are then indeed negligibly small. 



D. Evaluation and discussion of the rate 

We are now in the position to evaluate the rate-formula 
( pl| ) in terms of the master path x*p^{t). To this end 
we approximate in (| 9 2| ) and ( p3| ) the square brackets 
by 1, n eglect in ( |107D the term of order p*pj(tfc)^ and 
in ( |104| ) the last term (being also a correction of order 
By dropping the index of i/ we then can infer 
q) our central result for the instantaneous rate 



from 



r{t) = Vd aopt e-^-'/^ Koptit, D) [1 + OiD-')] (108) 



aopt := [4^r2 hm pl^.it)^ QlptitT^'^ (109) 



•-opt 



{t.D) 



k=0 



D 



X exp 



1 

D 



p*{t' + kTfdt' 



(110) 



The effect of our various approximations in deriving 
this result together with the corresponding "accuracy ex- 
ponent" 7 > in (108) will be discussed in Sect. IV. E 
below. Next, we analyze in more detail the pro pert ies of 
Kopt(i, D). By means of (|l]j6|,|7|) we rewrite ( pl0| ) as 



T^Opt 



K{t,ta) ^ 

Popt\^) ^ 



k=0 



D 

\2 rik 



X exp ■ 



2D 



C := 



-2 A„ r 



(111) 
(112) 



Since < C < 1 there is a competition in the sum (111) 
between the exponential terms which increase with k and 
the pre-exponential factors which decrease with k. One 
readily sees that the dominant contribution to the sum 
stems from a few fc- values around the real number fc, 
implicitly defined via 



p:p,it)^C''Iuit)^2D 



(113) 



Recalling that t stands here for tr and since neither 
Iu{t = tf) nor plpt{t = tf) (cf. ( [m] , ^ ) are small quanti- 
ties, it follows that k is, for small noise-strengths D, much 
larger than but, for sufficiently large tf — to, according 
to (|5^) also muc h smaller than_KYt = tf,to). Therefore 
the sum in (ITl) and thus in (IIC) can be extended to 
arbitrary integers k at the price of an error which is expo- 
nentially small in D, i.e., without actually affecting the 
accuracy exponent 7 in ( |108| ). As a further consequence 
of the fact that the dominant fc-values are much smaller 
than the upper limit K{t, tg) for large t — t^, we see that 
our formal approximation (pq) is indeed self-consistcntly 
satisfied. 

Next we notice that under the sum in 



110), the pre- 



exponential term is nothing else than the time-derivative 
of the expression in the exponential. By extending the 
sum over all integer fc-values as justified above we obtain 



1 

T 



t+T 



Koptit',D)dt' 



(fe+l)T 
t=kT 



= 1 — exp 



D 



Plptit'fdt' 



(114) 



Neglecting as usual errors exponentially small in D this 
leads us to the remarkable conclusion that 



1 /■*+^ 

^ J K,pt{t',D)dt' = 1 



(115) 



for all t and all (small) D. For t he time averaged rate 
(Q) we thus obtain from ( 10^ , 115| ) our central result p7| ] 



Dao 



pt' 



[1 + OiD-')] 



(116) 



It consists of an Arrhenius-type exponentially leading 
part with an "effective potential barrier" (/)opt and a non- 
trivial pre-exponential Z?-dependence. The two quanti- 
ties ckopt and (poot follow from the master-path x*p^(t) 
according to (|54| ) and ( 109 ). Thus they are independent 
of D but depend in a highly nontrivial way on various 
global properties of the deterministic force field F{x,t) 
in (^. In general, their explicit value can only be deter- 
mined numerically or by means of approximations. An 
exactly analytically solvable special case will be presented 
in Sect. V.A. 

We recall that for equilibrium systems, characterized 
by a time-independent force field F{x) = —V'{x) in (||), 
the escape rate exhibits an exponentially leading Arrhe- 
nius factor, which involves simply the barrier against es- 
cape of the static potential V{x), and a Z3-independent 
pre-exponential factor which depends only on local prop- 
erties of the potential at the barrier and the well jl ] , see 
also (161) below. The different structure of (116) is thus 
a consequence of the far from equilibrium situation cre- 
ated by the time-dependence of the deterministic force 
field F{x,t). 
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As announced in Sect. II. C, the time-averaged escape 
rate for the periodic force field F{x, t) can be identified 
with that of its supersymmetric partner force field jl^ ) 
for asymptotically weak noise D without any further re- 
strictions on F{x,t). The detaile d pro of of this highly 
non-trivial invariance property of ( 116 ) is carried out in 
Appendix D. 



Returning to the instantaneous rate (108), we see that 
i t ex hibits in comparison with the time averaged rate 
( 116 ) the additional time-dependent factor Kopt {t,D). 
The explicit evaluation of this factor requires the knowl- 
edge of one more global quantity, for instance of 



/3opt(i) hm t(*)e 



A„(t,t) 



(117) 



Note that due to relation ( |65|) the i-dependencyof this 
quantity is actu ally quite simple. Accor ding to (p^), this 
definition (117) allows us to rewrite (111) - with the 
range of k extended to arbitrary integers - as 



k— — OG 



D 



X exp 



2D 



(118) 



Besides l3opt{t) all other quantities in this expression 
are determined by local properties of the force field 
F(x u(t), t) alo ng the unstable periodic orbit. By exploit- 
ing (|6|j7|,|ll|) it follows that 



^opt 



it + T,D) 
it,CD) 



T-opt 



■^opt 
■^^opt 



{t,D) 
{t,D) 



(119) 
(120) 



Together with ( |113| ) and the obvious property < 
Kopt(t, D) < oo this completes our qualitative picture of 
the way in which T{t) oscillates around its average value 
f. 



E. The accuracy exponent 7 



In the following, we come to the determination of the 



accuracy exponent 7 in ( P^ , pl6| ). We wiU not elaborate 



here all the details of the rather involved calculations but 
restrict ourselves to the main steps. 

First of all, we recall that a contribution 0{D) is in- 
herited right away from formula (M). Next we have ap- 
proximated the square brackets in (|92|) by 1. For those 
/c-values which mainly contribute to the rate it can be 



inferred from ( |ll3| ) together with t = tj and ( p6[ ) that 



Plptitk^ = OiD) 



and hence with (106) that 



g:^dtk)Iuitk) = OiD) 



(121) 



(122) 



Since the integral in 
nitude as p*pt(tfc) from 
tal error we committed 




the same order of mag- 
we conclude that the to- 
is of the order 0{D'^), 



thus contributing on ce more a term of the order 0{D) 
in the rate- formulas (lOS , 116| ). The same conclusion can 
be drawn with respect to our approximating the square 
brackets by 1 in ( |93| ) and neglecting the OQ3*p^(ife)^) - 
term in (107) as well as the last term in (104). In other 
words, the relative error induced by all our so far made 
approximations is of the order 0{D). 

What remains is a closer inspection of the approxima- 
tion (|5^ ) for F{x,t) in the neighborhood of Xs,u{t)- One 
readily sees that actually only the approximation in the 
neighborhood of the unstable periodic orbit Xu{t) mat- 
ters in our quantitative evaluation of the rate; the basic 
reason for this is once more our assumption to —00 
in (p5|). In case that (^6|) happens to be an exact iden- 
tity in this neighborhood of Xuj t), t hen the total error 
committed in the rate-formulas (108,116) is thus of the 
order 0{D). Otherwise, a closer analysis of the relevant 
perturbative corrections shows that the error introduced 
via the approximation ( |5^ ) amounts to corrections of the 
order 0(j)*p^(tk)) in the rate-formula, i.e., of the order 
©(-v/D) according to (121). In other words, we can con- 
clude that 



^1 if F"{x^{t),t) = 
^ '1/2 otherwise. 



(123) 



In the case 7 — 1/2 it is important that in the global 
quantities (/iopt, Q^opt, /3opt(0 from (^,109,117) the long- 
time limits are made and the exact master path is utilized 
without any further approximations. If instead in these 
definitions any finite reference-time in combination with 
relations based on the approximation ( p6| ) were used, 
then this would introduce a possibly very small but nev- 
ertheless Z?- independent error an d so ruin the asymptot- 
ically exact predictions (108,116) in the weak noise limit 
L> 0. 

In cases for which (^6|) is not exactly satisfied in the 
neighborhood of the unstable periodic orbit Xu{t) and 
hence 7 = 1/2, it is in principle possible to calculate 
perturbatively the corresponding corrections such as to 
arrive again at a reduced relative error 0{D) in the so 
improved rate-formulas, though the actual calculations 
and the resulting expressions become very complicated. 
On the other hand, further reducing the 0(£))-error is 
even in principle rather problematic since it would require 
to go beyond the saddle point approximation in the path- 
integral approach from Sect. III. 

At this point it may also be worth to recall from 
Sect. IV. A that for any fixed (however small) D-value, 
the error 0{D) in (pl|), which is inherited by the final 
rate-formula, diverges as the amplitude of the time de- 
pendency of F{x, t) tends to zero, but also if its period T 
either tends to zero or to infinity. Thus, neither of these 
limits commutes with the limit D ^ 0. 



12 



F. The limits t ^ oo and D 



In the derivation of the rate-formula ( |108| ) we have as- 
sumed that all paths xl{t) which notably contribute in 
( pT| ) sojourn for a very long initial time-interval close to 
the stable perio dic orbit Xs{t), see (|85|). On the other 
hand, eq. ( 113 ) tells us that the amount of time which 
those dominant paths spend in the neighborhood of the 
unstable periodic orbit Xu (i) is roughly speaking of the 
order 0{\nl/ D). Both these conditions are compatible 
only if t — to substantially exceeds in order of magnitude 
\nl/D. In the physically relevant case, the noise-strength 
D is small but finite and this condition is well satis- 
fied after a comparatively short "tran sien t" time-period. 
Thus, strictly speaking, in (108) and (116), with decreas- 
ing _D-values also the lower limit of the admitted times 
t — to is tacitly assumed to slowly increase in proportion 
to Inl/D. 



We remark that our result (108) obvio usly remains pe- 
riodic in t for arbitrarily large t—to, see ( 119). T herefore, 
the restriction of the utilized basic formula (14) to values 
of t — to much s mall er than 1 /T does no longer apply to 
the final result ( |108| ), see also the discussion below (|l5|). 

In the physically less relevant case that t~tQ is kept at 
an arbitrary but fixed value and then D is made smaller 
and smaller, things become different. As pointed out in 
Sect. IV. A, for any finite initial and final times to and t = 
tf, there exists generically a unique absolute minimum 
xl^{t) of the action. For sufficiently small D the fco - 
term will thus completely dominate the sum in (|5l|), i.e., 

P*k (i/) e~'^''o(^"(*/)'*/)/-'^ 



nt = tf) 



[AnDQlitf]]^/^ 



(124) 



While most of the quantities on the right hand side of this 
result (including the index fco) still depend in a very com- 
plicated way on the time t = tf, no additional implicit 
-D-dependence is hidden. The most striking feature is the 
1/Vd pre-exponential behavior in comparison with the 
\/^ -scaling in ( p8|) ! 



Qualitatively, the crossover from (124) to (108), either 
as t increases or as D decreases, is clear: At some point 
the fc-dependence of the pre-exponential factors in ( ^l| ) is 
no longer negligible in comparison with the exponentially 
leading contributions and so the dominant /c-value moves 
away from fcg — ko{t,to) towards smaller values k k, 
cf. ( 113 ). At the same time, more that one term in the 
sum (|51) starts to notably contribute. 

Quantitatively, a leading order approximation follows 
along the same line of reasoning as in the derivation of 
( |108| ) from (|5l| ), except that in the approximation for 
the action (|91|), also contributions due to the deviations 
between x*f.{t) and its associated master-path a;*pt(i -I- 
kT) at times close to to have to be included, that is, 
the approximation ( ^ ) should be abandoned. The final 
result is again of the same form as in (1 
larger error than 0{D'^) and instead of (110) with 



but with a 



k=0 



X exp ■ 





" poo 






J. ^ 


J — oo_ 



D 



p:^,{t' + kTf dt' \ . (125) 



For moderate t — to oi extremely small D the exponen- 



tial in (125) depends very strongly on k and therefore the 
sum is dominated by a single term k — ko{t, to). Upon in- 
creasing t — to or D this strong /c-dependence of the expo- 
nential and hence the dominance of the fco-contribution is 
softened and the already discussed qualitative crossover- 
behavior is recovered. 



G. More general seeds xq 



So far, our rate-formula ( |108D is restricted to the case 
( ^ ) that the initial position xo at time to coincides 
with the stable periodic orbit Xs{to). As pointed out 
in Sect. II, one expects that for large enough times t — to 
the initial position xo should not matter, provided it is 
chosen inside the domain of attraction of a;s (t) . For suffi- 
ciently small noise-strengths D this is the case whenever 



Xo < Xu{to) 



(126) 



In the following, we analyze this intuitive expectation in 
some more quantitative detail. 

For arbitrary but fixed xo satisfying ( |126| ), the obser- 
vation from Sect. IV. A remains true, namely that only 
minimizing paths x^. {t) play a role in the rate ( pl| ) which 
closely follow a deterministic behavior i^. (t) ~ fX^k (^) ' ^) 
for most of the time. This requirement can be fulfilled 
in two basic ways and appropriate compromises thereof. 
The first possibility is that the path x^ (t) closely approx- 
imates a deterministic trajectory for a very long initial 
time- interval. During this time, a;^(t) approaches the pe- 
riodic attractor Xs(t) very closely and does practically 
not accumulate any action (36|3^). Consequently, one 
is back to the case (^2|) after an appropriate redefinition 
of the initial time to . Regarding the prefactor Qt (t) , 
one can, according to (^), approximate pl.{t) in (46) by 
zero. With the initial conditions ( ^7| ) one then recovers 
the same solution as in (^6|) except that in ( |63p^ ) the 
function As{t,to) is now defined as J^^ F'{xl{t'),t') dt' . 
Since x^(t) practically agrees with Xs{t) during a very 
long time interval, one sees that also with respect to the 
prefactor Ql{t) we are back to the case (|2|). As before, 
we may label such paths by low fc-values and their con- 
tributions to th e rate (E l j) are identical to those of the 
low fc-values in ( p^sHTToT 

The second possibility is that the minimizing path 
xl.{t) travels from its starting point xo immediately into 
the neighborhood of the unstable periodic orbit Xu(t) 
and then very closely follows this deterministic trajec- 
tory Xu{t) for the rest of its time. If xo is already close 
to Xu{to), such paths lead to a very small value of the 
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action in (36|3^) and thus will ultimately dominate the 
rate (|5l| ) if i — to is kept fixed and D becomes asymptot- 
ically small. This puzzling observation has lead to some 
amount of confusion in the recent literature . The 

resolution is that, much like in Sect. IV. F, things become 
very different for a small but fixed D in combination with 
larger and larger times t—to. The salient point is that the 
price to be paid for a long sojourn close to Xu{t) is a very 
smah p refa ctor pl{tf)/[Ql{tf)]^^'^ in ([5l|), as discussed 
below (107), namely of the order exp{— 2A„ • {t — to)}- 
As a consequence, the paths with low fc-values, as dis- 
cussed in the preceding paragraph, will dominate in spite 
of th eir unfavorable action. Ther efore , the rate-formula 
(108) appHes for any xq satisfying (|126D on condition that 



t-to > ct>opt/{2DXu) 



(127) 



This condition characterizes the asymptotic time regime 
for which the rate-formula (108) is valid in the case of 
a general initial condition. Even for rather small D, 
the preceding transient regime is typically confined to 
a few driving periods T, a s illustrated by the examples 
in Sect. V. Note that (127) comprises the condition from 
Sect. IV. F that t ~ to has to substantially exceed in or- 
der of magnitude In 1 /D. In other words, for a generic 
initial condition, cq. (127) is the only restriction for the 



rate- formula (108), apart from the exclusion of vanishing 
driving amplitudes and vanishing or diverging periods T. 



H. Summary from the practical viewpoint 

Given an arbitrary time-periodic force field F{x, t) that 
satisfies the condition in (H) , what are the necessary prac- 
tical (numerical or analytical) steps for an explicit quan- 
titative evaluation of the rate (108)? 

The first step is the determination of the stable and un- 
stable periodic orbits Xs{t) and Xu{t). An efficient way 
to do this is by evolving the deterministic dynamics for- 
ward and backward over a long time, respectively, with a 
reasonably well chosen initial condition. Once Xs^u{t) is 
known, the functions As^u{t,ti) from (63) and Is.u{t,ti) 
from (^l|,^3|) follow readily, with ti being an arbitrary 
reference time. 



The next step is the determination of the master path 
a;*pt(t). To this end, we chose an arbitrary but fixed 
time ts and a very small but finite positive number 
Axinin, characterizing the neighborhood of Xs{t) within 
which we are willing to accept the errors introduced by 
the approximation (p6|). We now consider the quantity 
Aa;*p(.(ts) as a parameter that may take values in the in- 



terval [Aa^mini 



Each such parameter value 



Ax*pj(is) yields a set of initial conditions 

Xlp,{t,) ^ Xs{ts) + ^xlp,{ts) (128) 

p:p,{ts) = /\x:p,%)/iSs) , (129) 

see ( |60| , p2[ ). With these initial conditions one then evolves 
x*pjf) and Popt(^) according to ( |4l]j4^ ). For a generic 



value of the parameter Aa;*pt(ts), the path a;*pt(t) will 
either reach Xu (t) after a finite time and then proceed to- 
wards a; = cxD or never reach x„ (t) and instead return into 
the vicinity of Xs(t) as t grows. By fine-tuning Ax*p^{ts) 
one has to find a path x*p^(t) which remains close to 
Xu{t) as long as possible, say until t = tmax- Upon vary- 
ing Ax*p^{ts) within [Aa;,„i„, Axmine^'''^ ^] the existence 
of at least one such path is guaranteed by the theory. A 
second solution, corresponding to a saddle point instead 
of a minimum of the action, is also to be expected (see 
Sect. V.A, below eq. ( |158| )). Further local extrema may 
coexist as well. Among them, the desired solution x*^^ (t) 
is the one with the smallest value of the action 



"opt 



2Is{ts) 



PUitfdt 



(130) 



see ( [43| , p4[|6l| , |76| ). By approximating t in (117) by t-^ 
we obtain 



/3opt(i)=P*t(imax)e^" 



(131) 



whence Kopt (i, D) from ( 111 81) follows with C from ( |ll2[ ) 
Finally, one chooses the initial conditions 



Q*{ts) - Isits)/^ , Q*,{ts) = is{ts)/2 



(132) 



see ( |99| , |100[ ), and then propagates Q%pt{t) according to 
( ^ ) until t — tmax, to obtain 

aopt = [4 7rr2p* (W)'Q:p t(imax)]~'/' , (133) 



13C 



131 



133) 



see (11091) . 

The accuracy of 0opt, /3opt(0i "opt from ( 
can be estimated by observing how little these quantities 
change if tmax is varied and if Ax^i^ is changed by a 
factor e'^"^. 

We finally note that the association of xlp^{t+kT) with 
the specific path xl{t) as in the preceding subsections 
does not play a role any more in the above described 
practical procedure. 



V. EXAMPLES 

In general, the explicit quantitative evaluation of 0opt , 
ttopt, and Kopt(t, D) in the rate- formula (108) is not pos- 
sible in closed analytical form. Exceptions arc piccewise 
parabolic potentials V{x) in conjunction with an additive 
sinusoidal driving (^). In Subsect. A, the simplest exam- 
ple with two parabolic pieces will be worked out and 
compared with accurate numerical results and with the 
approximation for small driving amplitudes from ||2^ . In 
Subsect. B we elaborate as a second example the case 
of a force field (||) deriving from a cubic potential V{x) 
along the lines of the numerical recipe from Sect. IV. H. 
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A. Piecewise parabolic potential 

We consider the force field from (^) witli a piecewise 
parabolic potential of the form 

A., 



V{x>Q) = ^[xl-{x-x^n 



(134) 



where x^ denotes the potential well (stable fixed point) 
and Xu the saddle (unstable fixed point), with the prop- 
erties 



Xs < , Xu> . 



The parameters 



A, < , A„ > 



(135) 



(136) 



characterize the piecewise constant curvatures and thus 
the time scales (Lyapunov exponents) of the determin- 
istic motion near the attractor Xs and the repeller Xm 
respectively. The force-field (||) then takes the explicit 
form 



F{x < 0,t) = \s(x -Xs) + A sinint) 
F{x > 0,t) = \u{x - Xu) + A sin(rit) 



(137) 



In particular, the quantities Xs^u in (134,137) are identical 
to those from (|68[ 



Requiring continuity at a; = we conclude from (137) 
that \sXs — \ii Xu- Selecting as independent model pa- 
rameters \s, A„, and the static potential barrier 



Ay := V[xu) ~ V{x,) (138) 
the fixed points be expressed through 



\2^V\\s\\u 
|As| + A„ 



(139) 



Turning to the determination of the stable and unsta- 
ble periodic orbits (Q) it is convenient to make a some- 
what stronger assumption than in (^ namely that both 
periodic orbits Xs^u{^") never cross the matching point 
x = of the two parabolic pieces of Vix), i.e. we re- 
quire that 



Xs{t) < < Xu{t) 



(140) 



for all times t. One finds that this property is granted if 
and only if the conditions 



A^ <(\lu + V?)xlu 



(141) 



are satisfied for both the 's'- and the 'u'-indices, and that 
the periodic orbits then take the explicit form 



A [As,„ sin(fi t) + cos(ri t)\ 



41) it follows that 



With the definitions ( |63| , 

^s,u{t, tl) = As,« • (< — ti) 
Is,u{t) = \Xs,u\ ^ ■ 



(142) 



(143) 
(144) 



Our next goal is the determination of the master path 
(t). To simplify the analytical calculations we restrict 



''opt 



ourselves to the case that the master path x^^^ (t) crosses 
the point x = exactly once, say at the time t = ti: 



t = u 



(145) 



The self-consistency of this assumption with the final so- 
lution fo r X* f (t) remains to be checked later. 

From ( |41p^ ) we see that both x*^^ (t) and p*pt (t) are 
still continuous a,t t = ti. For all other times t the re- 
lation ( |56| ) and the hence following conclusions are not 
approximations but exact identities since the force-field 
F{x,t) in (|137|) is by construction piecewise linear. By 
introducing ( |144| , ^5 ) and (|^) into (f|,|9|) we obtain by 
letting to oo and — s- oo for the master path the 
following two relations (one with index 's' and one with 
V): 



XsAtl) =Popt(*l)/As 



(146) 



These two equations for the two unknowns ti and p*p^(ti) 
imply with ( 142 ) the result 



tan(riti) = 

Popt(^l) = 



1 As A„, — ri^ 



n 



As + All 

AXs XuCoa{nti) 
r!(As + A„) 



(147) 
(148) 



We observe that the solutions ti of (147) are independent 
of A. Furthermore, there are obviously two solutions ti 
within every time-period T = 2Tr/n. We anticipate, that 
only one of them corresponds to a minimum of the action, 
and thus to the master path. Hence we fix ti uni quely 
(up to the usual degeneracy under t i—^ t + T) hy ( 147| ) 
in conjunction with 



AQcos{nti) 
As + All 



< 



(149) 



and show l ater, that this condition singles out the right 
solution of (147). [The case As -I- A ^ = ha s to be treated 
as limit As+A„ 0.] Combining (147- 149| ) it follows that 

Plp,{h) = A„x„ - \A\ |As| A„/i^2 > Q ^ (150) 
where we have introduced the definition 



(151) 



Note that XuXu in (150) may be rewritten in various 
equivalent forms according to ( |139| ) and that the last re- 
l ation p*t(ti ) > in (15C) follows as a consequence of 
(|T36|,pO,M6|). 
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Given ti and plp^{ti), the entire time evolution of 
the master path x*^f (t ) ca.n be readily inferred from 
(|^,|lH|j7|j7^,|7|) and (|l4|,|l4|) with the result 



Kpt(i) 



''opt 



(152) 
(153) 



where 's' is associated with times t < ti and 'u' with 
t > ti. All the general qualitative features discussed in 
Sect. IV. A. are nicely illustrated by this explicit example 
(|l5|,|53|). 

Finally , we hav e to check the self-consisten cy o f the 
solution (152,153) with our initial assumption (145), i.e., 
we have to verify that x*p^{t) is strictly positive for t > ti 
and negative for t < ti. In general, in doing so, a 
transcendental equation arises which has to be evalu- 
ated numerically. Without going into the details of the 
proof we further mention that one can show analytically 
that < As ,j is a sufficient but not necessary self- 
consistency criterion for ( |l45| ). On the other hand, it 
i s ob vious that the assumption Xs{t) < < Xu{t) in 
(140) is au tomatic ally cove red by the stronger require- 
ment (145). Thus, (140,141) is a necessary but not suffi- 
cient self-consistency criterion for (|l45 ). 

Introducing the above relations (152 |155 ) into ( p3| , p4[ ), 
we obtain for the action of the master path 



Oopt 



AV 



As A„(|As| -f A„) 



2AViy4 



1/2 



(154) 



For A — !■ or — > oo we thus rec over the static (un- 
driven) potential barrier AV from ( 138 ). The leading 
order corrections for small A decrease like |^| |2^. For 
any finite amplitude A and driving period T = 27r/f2 
the "effective potential barrier" 0opt is smaller than the 
static barrier AV and is monotonically decreasing both 
with increasing A and increasing T. Invoking the nec - 
essa ry b ut not sufficient self-consistency criterion (141) 
for (145), one can explicitly confirm that 0opt can never 
Ipl)) by dem onst rating that the 



become zero (see (^,| 

argument in the square-brackets in (154) is a lways pos- 
itive. If we had chosen the sol ution of (147) with the 
opposite inequality than in (149), then a plus instead of 
the mi nus sign in (152) would have been the consequence. 
Thus (149) is indeed the pertinent condition for singling 
out the solution which minimizes the action. 



By using (143) and (152) in the definition (117) of 
/3opt(i) we obtain 



/3opt(0=P:pt(ii)e-^"-(*-*^) . 

Turn ing to the prefactor Qopt('')i we see from 
(III) that 



1 



2IA., 



(155) 
j) and 

(156) 



for all tim es t < ti. Since F"{x,t) 
cording to (137), we can infer from (p 



= (A„ - Xs)6{x) ac- 
|) that the prefactor 



Q*p^{t) is continuous at t = ti while its derivative jumps 
from Qopt(^r) = to the value 



Qoptl^f ) 



|A.| 



''opt 



(157) 



where indicates the limit t — s- ti from above and 
from below. With these initial conditions, the solution 
of (p4) in the domain t > ti is straightforward, yielding 



(158) 



hm Q:pt(t)Kpt(t)' 



2 An 

Using ( 14^ , 15^ , 15^ ) one can show that p*pt(ti) - i;*pt(ti) 
i s id entical to the expression on t he left hand side of 
(149), so that (157) and thus ( |158| ) ar e pos itive quanti- 
ties. With the opposite inequality in (149) they would 
be negative, confirming once more that the latter case 
corresponds to a saddle point rather than a minimum of 
the action. Collecting e very thing, we are finally in the 
position to evaluate eq. (109) with the result 



^opt 



\A\ (r!2 + AsA„) + 



2AViA 



16 7r3|A|(^opt 



1/2 



(159) 



Once again, the fact that the argument in the square 
root is positive can be explicitly verified by exploiting 
t he n ecess ary but not sufficient self-consistency criterion 
(|T4l|) for (|l45|). 



B. Comparison of analytical and numerical results 

We have compared th e ab ove analytical predictions for 
the instantaneous rate ( |108[ ) with very accurate numeri- 
cal results in fig. 4. To this end, we have computed the so- 
lution of the Fokker-Planck-equation ( |l0| ) and then eval- 
uated the rate according to (|l3|), starting with an nar- 
row Gaussian initial distribution p{x,to) about the po- 
tential well Xs and then waiting until transients have 
died out, i.e. until T(t) has reached its T-periodic 
asymptotic behavior. In order to numerically evolve the 
one-dimensional time-dependent Fokker-Planck-equation 
( |To| ) one can employ standard parabolic partial differen- 
tial equation solving procedures in one spatial variable. 
We have adopted a Chebyshev collocation method ^] 
to reduce the problem to a coupled system of ordinary 
differential equations, which is then solved by standard 
numerical methods. By changing the various parameters 
of the numerical procedure, the typical relative errors of 
the numerical rates T{t) in our figures are estimated to 
be at most of the order of 10"'* for rates down to about 
10""'^°° and of the order of 10"'^ for rates down to about 

10^200^ 

The results in fig. 4 confirm for a representative set of 
parameter values that the agreement between the ana- 
lytical predictions and the practically exact numerical 
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results for the instantaneous rate r(t) indeed improves 
with decreasing noise-strength D. While the absolute 
values of T{t) and the location of the extrema strongly 
depend on D, the overall shape changes very little and 
does not develop singularities as D ^ 0. 

The corresponding time averaged rates (|ll6| ) are de- 
picted in fig. 5, exhibiting excellent agreement between 
theory and numerics even for relatively large D. The in- 
set of fig. 5 confi rms our analytical prediction that the 
relative error in (116) decreases asymptotically like D, 
see ( p^ . 

Finally, fig. 6. illustrates the dependence of the time 
averaged rate T upon the amplitude A of the periodic 
driving force. As expected, our theoretical prediction 
compares very well with the (numerically) exact rate, 
except for very small driving amplitudes A. The lat- 
ter discrepancy is in accordance with our discussion in 
Sect. IV.A and Sect. IV.E. 

We have furthermore included in fig. 6 a comparison 
with the analytical approximation for the time averaged 
rate f from Ref. p^ . By way of a matching procedure, 
involving the barrier region only, it is predicted p3| that 



r = rn 



2n 



(160) 



where Fg is the well-known Kramers-Smoluchowsky-rate 
in the absence of the periodic driving force B 



Fn = 



|y"(x,)F"(x„)|^/' _ 



AV/D 



27r 



(161) 



The leading order effect of an additive sinusoidal driving 
(^ , such that the associated periodic modulations of the 
potential barrier are small in comparison with the un- 
perturbed barrier AV, but not necessarily in comparison 
with the noise-strength Z?, are captured by the function 
s{(j)) in (161). It can be written for a general metastable 
potential V{x) in (|^) under the form ||2^ 



s{(j)) = yl (5 sin (/i — C cos ( 



S = S{xi) 
C = C{xi) 



dx sin ^57 
da; cos ( J7 



- dy 



X, V'{y) 
- dy 



, V'{y) 



(162) 



(163) 



with an arbitrary reference point xi G (xs,x„). Fig. 6 
confirms that this approximation from ||2^ is indeed com- 
plementary to ours in that it is very accurate for small 
driving amplitudes A but develops considerable devia- 
tions with increasing A. Those approximations have been 
omitted in figs. 4 and 5 since they are not valid in this pa- 
rameter regime and indeed are way off. 



C. Cubic potential 

As a second example we consider a force field (||) with 
a cubic metastable potential 



V{x) ^ --x'^ + - X, a,b>0 



(164) 



The stable and unstable fixed points Xg^u of this potential 
are given by 

= , x„ = - , (165) 
a 

with curvatures at those fixed points 

V"{xs)=h , V"(x^) = -h, (166) 
and a static potential barrier height 

&3 



AV ■.= V{x^)-V(x,) 



6a2 



(167) 



The time-dependent force field (^ takes the following 
form: 



F(x,t) = 



A sin(f7 i) 



(168) 



Since already the analytical evaluation of such a force 
field's periodic orbits is impossible, one has to recourse 
to numerical methods for the calculation of the quan- 
appt , and Kopt{t,D) appearing in the rate 



titles 



t'opt, 



expressions (108,116). A convenient numerical strat- 
egy for doing so has been discussed in detail already in 
Sect. IV. H. The s o obtained predictions for the time av- 
eraged rate (116) are compared in fig. 7 against precise 
numerical results for a representative set of parameter 
values. Note that these parameter values are quantita- 
tively very similar to those in fig. 6, hence also the rates 
as a function of the noise strength D are very similar. 
The agreement between the theoretical prediction and 
the practically exact numerical results is again excellent 
even for relatively large noise strengths D. However, in 
contrast to the piecewise parabolic case, the numerically 
accessible _D-values are still not smal l eno ugh in order to 
check the validity of our prediction (123) for the behav- 
iorof the relative error in the analytical approximation 

(HI). 



VI. CONCLUSIONS 

In our present work we have scrutinized by means of 
path integral methods the thermally activated escape of 
an overdamped Brownian particle over a periodically os- 
cillating potential barrier in the most challenging regime 
of weak thermal noise in combination with moderately 
strong and moderately fast driving. 

A first major result of our novel path-integral approach 
is the expression ( ^ for the instantaneous escape rate, 
which displays the suggestive general structure of "prob- 
ability at the separatrix times velocity" . The summation 
appearing in this expression reflects the fact that several 
local minima of the relevant action in the path-integral 
formulation of the escape problem notably contribute 
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to the rate. In contrast to the undriven escape prob- 
lem, giving rise to a (quasi-) Goldstone mode due to the 
(quasi-) time-translation symmetry, in our present case 
the paths corresponding to the local minima of the action 
are well separated and therefore admit a standard saddle 
point approximation of the path-integral for small noise 
strengths D. Pictorially speaking, by switching on the 
periodic driving, the continuous time-translation symme- 
try of the escape problem is reduced to a time-discrete 
one. 

Our present explorations indicate that from the practi- 
cal (technical) viewpoint, a path- integral approach which 
keeps an entire sum of possibly relevant contributions 
to the rate may be easier to handle than WKB-typc or 
quasipotential-type methods [^,^, which operate with 
the concept of a single exponentially dominating weak- 
noise contribution and a single pre-exponential factor, 
both of them typically of a non-analytic nature. 

The central result of our present paper represents the 
formula (108) for the instantaneous rate, supplemented 
by the result (123) for the "accuracy exponent" 7. The 
above discussed summation over the relevant local min- 
ima of the action resurfaces in all the equivalent alter- 



118) for Kopt{t,D) but drops 



native expressions ( 110| . Ill 
out (can be performed) in the time averaged rate (11£) 
due to the miraculous iden tity ( 115 ). 

The rate expressions (108,116) share the general 
Arrhenius-type structure of the exponentially leading 
weak-noise contribution with the typical form of an equi- 
librium (undriven) rate ( 161 ). However, both, the Ar- 
rhenius fa ctor and the pre-exponential contribution to 
the rates (|lO|,|ll|) now depend in a very complicated 
way on global features of the periodically oscillating po- 
tential (in contrast to the purely local properties AV = 
V{xu) — V{xs) and V"{xs,u) governing ( 161| )). Moreover, 
a non-trivial V^-dependence of the pre-exponential fac- 
tor on the noise strength D arises. 

For the time averaged rate ( |116| ) we have shown in Ap- 
pendix D that for asymptotically weak noise D an invari- 
ance property holds under the supersymmetry transfor- 
mation (Oq) without any further restrictions on the force 
field F(x,t). Such an invariance property can be estab- 
lished rigorously on very general grounds for force 
fields of the form F{x, t) = —V'{x) + y{t) and arbitrary 
noise-strengths D. 



The time averaged rate (116) displays a remarkable 
structural similarity with the rate-expressions obtained 
in | |50t for one-dimensional discrete-time systems in the 
presence of weak Gaussian white noise. While a general 
qualitative connection between these two different types 
of escape problems via some kind of stroboscopic map- 
ping is quite suggestive, the quantitative details are not 
so simple. Especially, the Gaussianity of the resulting 
noise after the stroboscopic mapping is crucial | |5l|] but is 
far from obvious for the rare but strong fluctuations 
(large deviations) which govern the escape events. 

The conditions for the validity of our rate-formulas 



are (127) and that for a fixed (small) noise strength D, 
extremely weak, fast, and slow periodic driving forces 
should be excluded. Especially, the weak noise limit D 
displays a rather intriguing non-interchangeability with 
the long-time limit < ^ 00 (see Sect. IV. F), and with the 
limits of asymptotically weak, fast, or slow driving. 

In general, an action-integral remains to be minimized 
numerically and an ordinary linear differential equation 
of second order for the prefactor to be solved (Sect. IV. H) 
before actual numbers can be extracted from our rate- 
formulas. However, for the special case of a sinusoidally 
driven, piecewise parabolic metastable potential this en- 
tire program can be executed in closed analytical from. 
This example retains all the typical features of more gen- 
eral setups and exhibits excellent agreement with high- 
precision numerical results (Sect. V). 

Conceptionally, our new path-integral approach should 
be of considerable interest for many related problems. 
Generalizations for higher dimensional systems and for 
non-periodic driving forces are currently under investi- 
gation [Q. Also the proper handling of the tantalizing 
weak, fast, and slow driving limits within a consistent 
path-integral formalism remains as an open problem for 
future research. 

Finally, the complicated dependence of the rate on the 
global details of the oscillating potential poses an chal- 
lenging inverse problem, namely to reconstruct the un- 
derlying force field from a given (e.g. measured) behavior 
of the escape rates. 
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APPENDIX A 



To prove the equivalence of (^9|) and ( |30| - |32| ) one 
first needs the Hessian of the discrete-time action 
Sn{xo, . . . ,xn) in (p4|), which is given by the (A^ — 1) x 
{N — l)-matrix 



n,m=l.2....,N-l 



( fli — ^1 

-61 02 



O.N-2 —bN-2 



(169) 



where 
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+ [ iF'{xl,tn)f + F«, i„) i„) ] Ai^ (170) 

&„ := l + ^^'(x;,i„)At . (171) 

For the prefactor Zjv(x*) in (^) the determinant of the 



matrix on the right-hand side of (169) has to be evalu- 
ated. This can be done by a standard procedure (cf. [^). 
The result is a linear two-step recursion relation for the 
principa l min or Q* , consisting of the first n columns and 
rows of (16£), of the form 



Q*n+l=anQ*n-bl_,Q:_, (172) 

for 2 < n < A^, with initial conditions 

Ql = l, Qi^ai. (173) 

Comparing (p9|j3^) with (|169|), we observe that — 



AtQ*pf, and due to the linearity of (172) we can conclude 
that Q* := AtQ"!^ als o obeys t hese equations. Therefore, 
using the definitions ( |l70|Jl7l| ), it is readily shown that 
Q* satisfies the required recursion relations (^0|) together 
with the initial conditions (|3l|). 

As a by-product, needed in Appendix B, we notice that 
by defining for 1 < n < A'^ 



IJ-n 



Q 



n+l 



Qn+1 



Qn 



(174) 



the linear two-step recursion relation ( |172| ) can be rewrit- 
ten as an equivalent non-linear one-step recursion 



Mn+l 



7 Ail 



ai 



(175) 



The actual quantity of interest then follows as 

N-l 



IN 



n 

n=l 



(176) 



APPENDIX B 



In the following we derive Eq. (pO[) by showing that 



2Qlitj) 



-F'ixf,tf) , 



(177) 



which gives together with ( |39|j4^ ) for t = tf the desired 
result. 

We work with the time-discrete version of the quanti- 
ties in (177) and consider in a first step the dependency 
of a minimizing path x* — x*{xf) on the end-point Xf 
for fixed to, tf, and xq. In order to complicate the no- 
tation as little as possible, we have left out the index k 
labeling the different paths . 



Since the initial point is kept fixed, we have that 
dx^/cLxf = 0. Further, we know from (p5| ) that for all 
n= 1,...,7V- 1 



dS 



N 



dXn 







(178) 



c*(2;/) 



for any x /-value, which implies after taking the derivative 
with respect to that 



dx*„, d^SN 



N-l 

y 

dXf dXjndXn 

m — 1 •> 



= 



Using (169), we thus get: 
dx ^ , dx 



bn- 



dx 



n+l 



dx f dx f dx f 

Introducing the new quantities ry„ by 

dxl 



Vn ■= bn 



dx* 



(179) 



(180) 



(181) 



(182) 



one obtains from (180) the one-step recursion relation 

bl 

Vn+l — fln+l 

The corresponding initial condition follows from ( |180D for 
n = 1 by taking into account the above mentioned fact 
that dx^/dxf — 0: 



ai . 



(183) 

Comparing (175) with (182 183| ) yields rjn — f^n for 



1 and therefore, using the definitions 
1 



(|171| , |174|181D , for n = AT 

[l + F'(x%_„tN-i)] 



IN 



dx 



N 



Q 



N-l 



dx 



(184) 



N-l 



In a next step an exphcit expression for dx*^^ / dx*^ 
in terms of well-known quantities has to be found. This 
can be achieved by taking the second derivative of the 
discrete-time action S'jv(x*(a;/)) of the same minimizing 
path X* as above with respect to the end-point x f . With 
(178,179), and the boundary conditions (26) we find 



d^SNi^*{xf)) 

dx} 



N 

E 



dx* 



'AT 



dx*j^ dxndxN 



n=l "-^N 

and using (|2j) we can conclude that 



c*(a;/) 



dx 



N 



dx 



N-l 



l + AtF'{x*j^_^,tN^i) 
1-2 At £^SNi^*{xf)) 



Together with (|184| ) we thus arrive at 



_ [l + AtF'{x*^_„tM-i)y 



l-2At^SN{^*{xf)) 



(185) 



(186) 



(187) 



which with ( p3| ) yield s in the continuous-time limit the 
searched for relation (177). 
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APPENDIX C 



In this appendix we derive ( |107D for < G [tu,tf], where 
the order of magnitude refers to the asymptotics with 
respect to k. 

As discussed below (83), the differences 



Sxlit) :=4(0-.<pt(i + fcr) 

Split) -.^pm-plptit + kT) 



(188) 
(189) 



rapidly decrease with increasing index k uniformly on 
the entire time- interval [io,^/]- Our first conclusion is 
that the time-point i„ at which a;*(t) enters the neigh- 
borhood of Xu{t) depends itself on the index fc, basically 
it decreases like —kT, see fig. 3. On the other hand, the 
distance Aa;^(i„) at which the path enters this neighbor- 
hood is, by definition, fc- independent. The corresponding 
momentum p^(iu) is not strictly fc-independent, but ap- 
proaches an asymptotic fc-independence for large k. Sim- 
ilar conclusions apply for the time tg at which j:* (t) leaves 
the neighborhood of the stable periodic orbit. 
Next we can conclude from (^j76) that 



Sxlit) = Split)I,it)~plit)e^^^<~'''o)l,it„) . (190) 

Within the approximation ( |85| ) it follows that Splits) — 
Sxlits)/Isits)- With these initial conditions at t = tg, 
the small perturbations SxHt) and Split) are then prop- 
agated according to ( p| , ^2| ) until t — tu- In linear order of 
these small perturbations it follows that SpHtu) / SxHty) 
is an asymptotically fc-independent constant, which, how- 
ever, depends on all the details of the force field F(x,t) 
along the crossover se gmen t of the master path x*p^(t). 

The counterpart of (|l90| ) in the neighborhood of a;„ it) 
follows along the same line of reasoning, reading 

Sxlit) = ^Split)Ut) +p^(t)e-2A.(t/,t) j^(i^) . (191) 

Replacing on the right hand side e~^^=(*-f by Po^titf + 
kT)"^ /plp^it + kT)"^ acc ordi ng to ( |6l| , p5| ), choosing t ~ i„, 
and making use of (p6,189) we can infer that 



Sxlitu) + Splitu) luitn 



Pl{tu)Iuitf) 
\Plitu) - Splits)]- 



(192) 



As we have just pointed out, the quantity (5x^(t„) 
is proportional to SpHtu) with an asymptotically k- 
independent proportionality constant that depends on 
the details of the force field Fix, t) along the crossover 
segment of the master path x^^^it). In the generic case, 
this proportionality constant is thus not expected to coin- 
cide with —luitu) since the latter depends on the behav- 
ior of Fix, t) along the unstable periodic orbit x„(i) only. 
Conseq uentl y, both SxHt^) and SpHtu) on the left hand 
side of ( [19 2D are, with respect to their fc-dependence, of 
the same order of magnitude as the right hand side. Since 



plitu) is asymptotically fc-i ndep endent and Split u) tends 
to zero, we can infer from (|l92| ) that 



Sxlitu) = Oip*itk)') , Splitu) = Oip*itk)') 



(193) 



With the initial conditions (^,100) it follows from (|4£ 
that the relative difference between QHtu) and Qoptitu- 
kT) sc ales as a function of k like SxHt^) and Split u 
With (|93l) this implies that 



Qlitu) = Qlptitu + fcT) [1 + Oiplp,itk)^)\ . (194) 

A similar relation follows for QHtu) and thus for glitu) 
(cf. (^)), namely 

9l{tu) = glp.itu + fcT) [1 + Oiplp.itu)^)] . (195) 



Finally we conclude from (101) that 

Ql(t) 



^ opt 



Qlitu) 

it + kT) Q*itu + kT) 



1 - glitu) lujtu.t) 
l~g*,ptitu + kT)Iuitu,t) 



(196) 



Like for Ax^(t„) and plitu) (see below ( [189D ) one can 
convince oneself tha t also g titu) is asy mpto tically fc- 
independent. With ( |l9^Jl95| ) the result ( |l07|) then fol- 
lows from dl9q). 



APPENDIX D 

The purpose of this appendix is to verify that our ex- 
pression ( |11(^ ) for the time-averaged rate is invariant with 
respect to the supersymmetry transformation (|l6|). To 
this end, we first note that the path defined via 



^opt{ t) 

Plpti-t) 



(197) 
(198) 



^opt (t) 

Plptit) 

satisfies the Hamilton equations ( ^ , ^2| ) for the super- 
symmetric partner field Fix,t) from eq. (p^. Since 
the periodic orbits of this new force field are given by 
Xsit) = —Xui~t) and Xuit) — —Xgi —t) (see Sect. I.C) 
one can readily see that i*pj(t) from (197) also obeys the 
boundary conditions ( |5^ ) in the relevant limit t/ — to ~* 
CO. Hence we have found (up to the usual degeneracy 
with respect to time shifts by arbitrary multiples of T) 
the unique solution of the supersymmetric partner vari- 
ational problem ( H ). Inserting p*pj(t) from ( |198| ) into 
the definitions ([43|) and (^J) then leads to the following 
result: 



(199) 



Somewhat more elaborate considerations are necessary 
in order to establish a corresponding identity for the prcf- 
actor ckopt in (116). To this end we first consider two arbi- 
trary but linear independent solutions Qiit) (i — 1, 2) of 
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the prefactor equation ( }46D for Q*^ (t) . One can then eas- 
ily verify that the prefactor Q*t jt) , which moreover has 
to fulfill the initial conditions in the limit to — > — oo, 
is given by 

n* r QiW Q2{to)-Qi{to)Q2it) 



with the Wronskian 

W{t) ■.= Qi{t)Q2{t)-Qi{t)Q2{t) 



(201) 



Due to ( [iq ) one can infer that 

W{t) = 2W{t)F'{xl^,{t),t) . (202) 

With help of the Hamilton equation (|4|) it follows that 

pl^,{tfW{t) =con^i. (203) 

Turning now to the supersymmetric partner problem, 
it is readily seen that one obtains via Qi{t) :— Qi{~t) 
two linear independent solutions of the prefactor equa- 
tion ( ^ ) for the su persy mmetric partner field ([l6 ) and 
the path given by (197). Thus we can use ( |198 ) and 
(200,201) (with tildes) to establish the identity 



t—*oo ^ '■ 

X [Qi{-t) Q2hta)- 



-W{-h) 
i{-h)Q2{-t)] . 



lim 



(204) 



According to ( |203| ) we can now 
rewrite pl^,{-tf /W{-to) as pl^,{-tof /W{-t). Re- 
placing t ~f — tp and vice versa one can then conclude 



with help of (g00|) that 

}}^pl^t(tf Qlptit) = }}^P*opt{t)' Qlptit) ■ (205) 
Hence we finally obtain 

ttopt = fiopt • (206) 

This, in combination with ( [199D, p roves our proposition 
that the time-averaged rate ( |116| ) is invariant with re- 
spect to the supersymmetric transformation (p^). 
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FIG. 2. The supersymmetric partner potential 

V{x) := —V{—x) of the potential V{x) from fig.l. 
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FIG. 3. Solid: The paths xlii), k = 0, . . . , K(t;, to) = 3 
which minimize the action (MJs;^) with boundary conditions 
(^sj). Dashed: The associated "master paths" Xopt (t + fc T) , 
implicitly defined via (^^. Dotted: Stable and unstable peri- 
odic orbits Xs{t) and Xuit) from (^). In this plot, tf — to has 
been chosen rather small. Astf — to increases, more and more 
intermediate paths x1{t) appear which better and better agree 
with their associated master paths x*p^{t-\-kT). The depicted 
curves have been obtained for the additively driven piecewise 
parabolic potential (^ 134) with parameters A = 0.5, Q = 1, 
As = —1, \u = 1, AV = 1, to = —12, tf = 7.5 (dimensionless 
units). 



FIG. 1. Sketch of a typical metastable potenti al V (x) in 
eq.(^. Plotted is the piecewise parabolic example (134) with 
parameters AV — 0.9, As = —0.6, and Au = 0.3 in arbitrary, 
dimensionless units. 
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FIG. 4. Instantaneo us ra te T{t) versus time t for the piece- 
wise linear force field (137) with parameters — Xg = —1, 
Xu = Xu = 1, ~ 1, and A — 0.5, corr espondi ng to a 
static {A = 0) potential barrier AV = 1 in (|l38|,E39|). Solid 



hue: Analytical prediction ( |l08| , [l 1 ^ll 5C| , |l 5 lyi 54[1 S^Jl 59[ ) 
by neglecting the 0{D"')-term in ( p_08| ). Dashed line: 
High-precision numerical results, obtained as described in 
Sect. V.B. 
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FIG. 6. Time-averaged rate F versus driving amplitude 
A for the piecewise linear force field (137) with parameters 
Xs — \s = —1, Xu = Xu = 1, O, — 1, an d D — 0.05. Solid 
line: Analytical pre diction (116,154,159) by neglecting the 
0(D" '')-te r m in (tllq). Dotted line: Theoretical approxima- 



tion (16C)-(163) according to Ref. [23] Crosses: High-precision 
numerical results, obtained as described in Sect. V.B. Inset: 
Magnification of the small- ^4 regime. 
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FIG. 5. (a) Arrhenius plot of th e time-averaged rate f for 
the piecewise linear force field (137) with the s ame parame ters 
as in fig. 4. Solid line: 
neglecting the O(D^) 



Analytic al pr ediction (116 154 159) by 
-term in (116). Crosses: High- precision 



numerical results, obtained as described in Sect. V.B. (b) Rel- 
ative difference between analytical (F) and numerical (Fnum) 
rate. 




FIG. 7. Arrheni us p lot of the time-averaged rate F for the 
cubic potential (^,164) with parameters a = l/%/6, b = 1, 
A = 0.5, and Q = 1, corr espo nding to a static {A — 0) poten- 
tial barrier AV = 1 in (167) and curvatures — 1 
in (166). Solid line: Analytical prediction from (116) with- 
out the C'(D''')-term by adopting the calculational procedure 
from Sect. IV. H. Crosses: High-precision numerical results, 
obtained as described in Sect. V.B. 
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